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Abstract. We consider the distribution of the value of the opti- 
mal fc-assignment in an m x n-matrix, where the entries are inde- 
pendent exponential random variables with arbitrary rates. 

We give closed formulas for both the Laplace transform of this 
random variable and for its expected value under the condition 
that there is a zero-cost k — 1-assignment. 



1. Introduction 

Let M be an m x n matrix. A fc- assignment is a collection of k entries 
in the matrix such that no two are in the same row or the same column. 
The value of a fc-assignment is the sum of its entries. A fc-assignment 
is called optimal if its value is no larger than the value of any other 
fc-assignment. If the entries of the matrix M are random variables then 
so is the value of the optimal fc-assignment, here denoted minfc(M). 

The study of the optimal fc-assignment has been pursued by re- 
searchers from different fields and with different random variables as 
entries in M. The main focus has been to estimate the size of the 
expected value of minfc(M). For references and more details on the 
history see [HS98] or [TWHSj . 

In 1998 Giorgio Parisi |P98j conjectured that if M is a fc x fc matrix 
with independent exponential random variables with rate 1, exp(l), 
then the expected value of the optimal fc-assignment is 

E[min fc (M)]=l + i + i + ... + l 

Two very different proofs of this conjecture were announced simul- 
taneously in March 2003, |1AV03| |NPS03j . The beautiful conjecture of 
Parisi inspired much work on exact formulas and many different gen- 
eralizations where studied [AT^OH MJBM IUS981 IUS021 lEESOl . LW00 , 



IN02j . In |LW03j a formula for the expected value is given when the 
matrix entries are exp(l) or 0. 

In this note we investigate the problem from a different extreme. 
We allow the rates of the exponential random variables to be arbitrary 
positive numbers. We include the infinity as a possible rate, which 
corresponds to the entry being constant zero. We prove exact formulas 
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for E[mirifc(M)] and for the Laplace transform L(min^(M)) under the 
assumption that E[minj( c _i(M)] = 0. Formulas for completing a k — 1- 
assignment were considered previously in jCS02j in the case when all 
rates are equal to 1. The Laplace transform for some special cases 
when the rates are all equal to 1, have been determined in (3-S02 . 

We will first prove the slightly easier result on expected value and 
then compute the Laplace transform for the entire distribution. In 
theory one should be able to deduce the first result from the second 
but it seems easier to compute them separately. 

2. Preliminaries 

As in jLWOOl ILW03| the concept of row and column covers of zeros 
will be important. The set of zeros will be denoted Z . We will consider 
sets of rows and columns in the m x n-matrix. A set a of rows and 
columns is said to cover Z if every entry in Z is on either a row or on 
a column in a. A cover with k — 1 rows and columns will be called a 
k — 1-cover. For many readers it might be convenient to translate the 
matrix to a bipartite graph. In that setting our covers are so called 
vertex covers. 

Given a set of rows and columns a let the rectangle R(a) be the 
part of the matrix not covered by a. If a is a k — 1-cover of the zeros 
then the corresponding rectangle R(a) will be called critical. Let Qk,z 
be the set of all critical rectangles in M. We define a partial ordering 
on Q kyZ by letting, for R(a) 1 R((3) G Qk,z, R(ot) < R(/3) if the set of 
columns in a is a subset of the set of columns in /3, and the set of rows 
in a is a superset of the set of rows in (3. 

Recall that a random variable X ~ exp(a) if P(X > t) = e~ ta . Then 
a is called the rate and E[X] = 1/a. If Xi, . . . , X n are independent and 
Xi ~ exp(a!j) then E[minjXj] = Let the rate of a rectangle 

R be the sum of the rates of the individual entries. We denote it I(R) 
and note that < I(R) < oo. 

From matching theory we have the following lemma that will be 
important. 

Lemma 2.1. Assume that Z contains k — 1 zeros in independent po- 
sition. Then Qk,z is a lattice. In particular it has unique maximal and 
minimal elements. 

From Theorem 2.9 in |LW00j we cite the following useful fact. 

Lemma 2.2. Suppose that a row r belongs to every k — 1-cover of the 
zeros Z in the matrix M . Then every optimal k-assignment contains a 
zero from row r. This means that we can remove row r from M and ob- 
tain a matrix M r with the property that E[min^(M)] = E[minfc_ 1 (M r )]. 

Recall that the incidence algebra over a poset is the algebra of func- 
tions defined on the intervals in the poset, see [H]. The product in the 
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incidence algebra is defined by convolution. One way to visualize the 
incidence algebra over a poset Q is to let a function / be represented 
by a matrix F with rows and columns indexed by the elements in Q. 
The element on position (a,/3) in F is f(a,/3). The rows and columns 
must be arranged according to a linear extension of Q. Every such 
matrix is upper triangular and multiplication of functions corresponds 
to multiplication of matrices. 

3. The expected value 

Define a function in the incidence algebra of Qk,z by f(R(a), R(j3)) = 
I(R(a) fl R(/3))- Since I(R(a)) > for all critical rectangles, / is 
invertible in the incidence algebra. 

Theorem 3.1. Let M be an m x n matrix with all entries being either 
zero or independent exponential random variables with arbitrary posi- 
tive rates. Assume that the set of zeros Z in M contains k — 1 zeros 
in independent position, i.e. E[min/ C _ 1 (M)] = 0. 
Then 



E[min fc (M)]= £ f-\R(a),R(j3)), 

R(a)<R(/3) 

where the sum is over all intervals in Qk,z and f" 1 is the inverse of f 
in the incidence algebra. Equivalently we can write 

E[min fc (M)] = 

, ( s-i) i(Ri n gjJCgg n R 3 ) ■ ■ -jjRs-i n R s ) 



I(RMR 2 )--I(R S ) 

where the sum is taken over all non-empty chains in Qk,z- 

Remark 3.2. The second formula runs over all chains in Qk,z which in 
the worst case has size of order k\. The first formula is a computational 
and conceptual improvement for large k. It involves taking the inverse 
of a matrix indexed by the elements of Qk,z which is exponential in k 
in the worst case. 

Proof. The equivalence follows from Lemma f5. II below. The proof will 
be by induction over k. The theorem is certainly true for k = 1. With- 
out loss of generality we may assume that the entries (1, 1), (2, 2), . . . {k— 
1, k — 1) are zero entries. We may also assume that the maximal rec- 
tangle -R(7) in Qk,z corresponds to the cover 7 consisting of columns 
1, . . . , k — 1. If this is not the case, then there is a row i that belongs 
to every k — 1-cover. This implies by Lemma (2.21 that E[min/ C (M)] = 
E[minfc_i(M i )], where M l is obtained from M by completely removing 
row i. Since these matrices have the same Qk,z the result is clear by 
induction. 
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Note that any chain in Qk,z not ending with R(j) can be augmented 
with i?(7) at the end. 

We now use the same recursion procedure as in |AS02j and |LW00j 
corresponding to the optimal cover 7 consisting of the first k — 1 
columns. That is, let X be the minimum of all the exponential ran- 
dom variables in #(7). Note that exactly one entry in #(7) will belong 
to an optimal ^-assignment and that E[X] = l//(i?(7)). Subtract X 
from all entries in R("f) and a new zero will occur at the position of the 
minimum. All other entries will be unchanged in distribution by the 
forgetfulness of the exponential distribution. Let My be the matrix 
obtained when position in M has been replaced with a zero. 

Let Ki be the intersection of R{^) and row i. If the minimum is in 
Ki, i > k — 1, then we get a zero fc-assignment in My. If the minimum 
is in Ki, i = 1, . . . , k — 1, then row % has to be in every optimal cover 
of Mij and as above we can remove row i and this case is done by 
induction. Let again M l denote the matrix with row % removed from 
M. Also let R l be the column maximal cover of M\ The poset of 
k — 1-covers of M % will be the induced subposet of Qk,z of elements 
< R\ 

The probability that the zero occurs in Ki is 77^^) an d we get 

E[,„ lM M)] = + g ^K-.^)]. 

Which by induction becomes 



'"' l(K,) ^ , /(ft nR 2 ) ■■■l(R,- l nfl. 



Change the order of summation to get 



1 



I(Rj) 



v (8 _ 1) /(ig 1 nig 2 )---/(fl s _ 1 niQ v 7(2^ 



i(r 1 )---i(r 8 ) ^ i(R(i)y 

Rl<-<R s <R(-y) K ' V J i(£rowset(R a ) v KUJ 

Here rowset(R) denotes the set of rows of M that intersect the rectan- 
gle R. Now 

v I{K) I(R(i)\R s ) _ I(R S n fl( 7 )) 



since the entries are independent exponential random variables. The 
theorem follows. □ 
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4. The Laplace transform 

We can in fact use the same proof technique to get the stronger result 
determining the Laplace transform of the distribution of min^M). 

Recall that the Laplace transform of a random variable X is L(X, t) = 
Fj[e~ tx }. It has the following well-known properties. 



L(X + Y,t) = L(X, t)L(Y, t), if X and Y are independent 

Given random variables X±, . . . ,X S and probabilities p±, . . . ,p s , de- 
fine the random variable I to take value % with probability pi, indepen- 
dent of Xi, . . . , X s . Then 

s 

L(X I ,t) = ^ /Pi L(X i ,t). 

i=l 

In this situation we will need the special case L(0, t) = 1. 
For a critical rectangle R we will use the notation 

(j)(R,t) = L(min 1 ( J R)) - /(/?) 



I(R) + t 

As for Theorem Kill we give two statements of the same formula using 
Lemma f5. II Remark 13.21 applies also here. 

Theorem 4.1. Let M be an m x n matrix with all entries being either 
zero or exponential independent random variables with arbitrary posi- 
tive rates. Assume that the set of zeros Z in M contains k — 1 zeros 
in independent position, i.e. E[min/ C _ 1 (M)] = 0. 

Then we can write L(minfc(M)) = ^2 RG Q k c R (M)<p(R,t) . Further- 
more 

c R (M) = a R ■ b R , where 

{ _ 1)s j^l(RinRi + i)-I(R) 



I(Ri) - I(R) 

Ri<R 2 <-<R s =R i=l v ' K 1 



and 



" = »*<S<...<». ( ' .5. m=m ■ 

where the sums are taken over all chains containing R in Qk,z- 
In the generic case rewrite this as 



c r {m)=[ 9r(r(<*),r))'[ E 9 R \R,Rm 
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where the sums are over elements in Qk,z an d where g 1 is the inverse 
in the incidence algebra of 

[l, if Ri = Rj = R 

g R {R i ,R j )= h, ifRi^Rj 
[I{Ri n Rj) - I(R), otherwise. 

Proof. The proof is by induction over k. We use the same notations as 
in the proof of Theorem 13.11 and compute the Laplace transform using 
the same recursive step. 

L(mm k (M),t) = 4>(R(j),t) + E jf^) " 

where p is the probability that the minimum is located so a zero cost 
/c-assignment occurs. We now decompose this product by the method 
of partial fractions with respect to t. This proves the first claim and the 
uniqueness of the coefficients c R (M). The function g R is invertible if 
I(R) ^ I(R(a)) for all R(a) 6 Q KZ such that R < R(a) or R(a) < R. 
The equivalence of the two formulas for the coefficients then follows 
from Lemma f5. II 

After the decomposition by partial fractions follows an extraction of 
the terms involving <j)(R, t) which gives 

Assume that R ^ R(j) and we may inductively write 

^ I(Kj) 
t J(i2( 7 )) -I(R)' 



1 J -LA I(Ri)-I(R) ' 

R=R 3 <Rs+i<-<Ru i=s+l v J/ v 7 

where i? n does not intersect the row z. Change the order of summation 
to get 

c R {M)=a R {M) (-l) ( - s+1) - 

R=R s <Rs+i<-<Ru 

A )(Ji,n6 H )-;(«) v /(A-,) 

and proceed as in the proof of Theorem 13.11 Here colset(R) denotes 
the set of columns of M that intersect the rectangle R. Note that 
trivially a R (M) = a R (M l ). If R = R(j) the recursive step uses the 
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row maximal rectangle instead, that is the smallest element in Qk,z to 
compute a#(M). □ 

5. A LEMMA 

Lemma 5.1. Let f be a function in the incidence algebra over a poset 
P and let a < (3 G P be arbitrary elements. The inverse of f can be 
written as 

f -l (n o\ V" ( 1 ^ 7l ' 72 )^ 72 ' 73 ) " " " [(jfzll 7 *) 

J 1 ,P) < ^ < 3 [ } /(7i,7i)/(72,7 2 )---/(7 s ,7 s ) ' 

a=7l<72<---<7s=P 

where the sum is over all chains in the interval [a, ft] beginning in a 
and ending in (3. 

Proof. Let F be the upper triangular matrix corresponding to / as 
described in the preliminaries. Let D be the diagonal matrix with 
the values of 7(7,7) on the diagonal and zeroes elsewhere. Let N be 
the nilpotent matrix that agrees with F at all positions except on the 
diagonal, where N has zeros. We can then write F = D + N and we 
can easily verify that 

F- 1 = IT 1 - D^ND- 1 + D^ND^ND- 1 

The matrix D~ l N will have zeros on and below the diagonal. Thus the 
sum is finite. Since inverting / in the incidence algebra is the same as 
inverting the matrix F the lemma follows. □ 

Remark 5.2. The main theorem in |LW03j when the exponential ran- 
dom variables all have rate 1 or infinity has a reformulation in jLWOOj 
in terms of the Mobius function of a certain poset called P. That poset 
is different and much larger than Qk,z- To be more precise the atoms 
in P are the elements in Qk,z- All our efforts to join the two theorems 
to one for completely arbitrary rates have so far been fruitless. 
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